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ON THE VARIATION OF THE ARBITRARY AND GIVEN CONSTANTS 
IN DYNAMICAL EQUATIONS* 

BY 

ERNEST WILLIAM BROWN 

Introduction. 

The method of the variation of arbitrary constants is a well-known device for 
obtaining the solution of a set of differential equations when we know that of 
another set which differs from the former in a suitable manner. Especially is 
it useful in most of the problems of celestial mechanics where methods of approxi- 
mation have to be adopted. The first approximation once obtained, it gives 
the second approximation with ease and certainty, and the same idea can be fol- 
lowed out for the higher approximations if further accuracy is desired. 

In the class of problems under consideration here, the equations of motion 
can be put into the form 

d 2 x. BF 
df ~ Bx. 

or into the form, known as canonical, 



(i=l, 2, ••■«), 



dx. BF dX. BF 

« J. 

dt ~ BX. ' dt ~ 3x~ 



(»•=!, 2, 



where F is a function either of the x i or of the x i , X. , respectively, and may 
or may not contain the time explicitly. Suppose that the problem has been 
solved so that the x i , X i are expressed in terms of 2n arbitrary constants and 
the time, and that it be required to find the solution of another problem in which 
F + H replaces F. The new problem can be stated by supposing the arbitrary 
constants in the first problem to be variable and finding the differential equations 
satisfied by them. 

The function F will usually contain, besides x., t, a number of given con- 
stants which will naturally be present in the solution, and it may happen that 
some of these constants are really variables which, in order to get a first approxi- 
mation, we have treated as constants. It is, of course, not possible in general 

* Presented to the Society uuder a different title April 25, 1903. Received for publication 
April 20, 1903. 
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to solve the equations as though they were constant and then substitute their 
variable values in the solution. If we follow the method mentioned we may ask 
what variable values are to be given to the arbitraries in order that the variable 
values of the given constants may be substituted in the solution. Thus instead 
of adding R to F the new problem is that obtained by replacing some of the 
constants present in F by functions of the time ; it is chiefly this problem that 
I propose to consider in the following pages. 

In order to do so it is necessary to consider in some detail the relations satis- 
fied by the derivatives of the dependent variables, not only with respect to the 
arbitraries but also with respect to the given constants. In the course of this 
investigation, generalizations of certain theorems which I gave in 1897 * (which 
included generalizations of Adams' theorems connecting the mean motions of 
the perigee and node with the constant term in the expression for the moon's 
parallax) are obtained and some new relations between the various constants are 
added (§§ 1-4). In § 5 the equations for the variable values of the arbitraries 
corresponding to variable values for the given constants are obtained. It there 
appears that Newcomb's theorem (see art. 24) is a particular case of a general 
theorem which applies to all equations of dynamics satisfying the conditions 
laid down. In § 6 applications of the results obtained are made to the lunar 
theory and in particular to the effects which variable values for the given con- 
stants produce on the motions of the arguments. 

§ 1. The Equations of Variations. 
1. Let (a;, y, z) be the coordinates, referred to fixed axes, of a particle mov- 
ing under a force-function F. The equations of motion are 

(1) 

The solution of these equations involves a certain number of constants : of these, 
six are arbitrary and the others are known and present in F. 

Denote by f, g any two of these Constants and differentiate (1) with respect 
to /, g successively. We find 

d^dx_d 2 ^dx d 2 F dy d 2 F dz d^BF 
dfdf~~dx Y df + ~Bx~By df + Bx^Bz df + df dx ' 

d 2 dx d 2 Fdx d 2 F dy d 2 F dz d'dF 

. ^ I *L _J_ L 

dt 2 dg dx 2 dg dxdy dg dxdz dg dg dx ' 
with similar equations for the other coordinates; the operators d'/df, d' /dg 



d 2 x BF 


d 2 y dF 


d 2 z dF 


dt 2 ~ dx ' 


dt 2 ~ dy ' 


df ~ dz 



* On certain properties of the mean motions, etc., Proceedings of the London Mathe- 
matical Society, vol. 28 (1896), pp. 143-155. 
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denote derivatives with respect to f, g only in so far as they occur explicitly 
in F. 

Multiply these equations by dx/dg, — dxjdf respectively, and add the cor- 
responding results for the other two coordinates. Denoting by S a sum of simi- 
lar expressions with respect to the three coordinates, we obtain an equation which 
may be written, 

d fdx d' 8F dx d' 8 FY d d'F d d' F 

( 2) dtV' g > = ^{dg df~dx~ df dg ~8x~ ) = dg ~df~ ~ ~df~dg' 
where 

fdxd'dx dxddx\ 

(3) (f^)= S {dgdtdf-drdtdg)--^^- 

If h be another of the constants we have, by direct differentiation, 

(4) ^(*.A) + |(*,/) + £(A*)-0. 

2. The constants will now be separated into two classes. Suppose that the 
solution is such that F and each coordinate is expressible by means of sums of 
cosines or sines of sums of multiples of the angles 

b t t + p i = te 4 , 

where h i and the coefficients of the periodic terms depend only on one set of 
constants a. , the other set being /3 f , and the a i , /3 { being all independent of 
one another ; these will be respectively called the a-set and the /8-set. It is sup- 
posed that F is given in terms of as, y, z and the known constants, involving 
the latter in the form just given. It is also supposed that the expressions for 
F and for its derivatives contain no term independent of x,y, z. Also the 
constant of length is to be such a function of the arbitrary constants of the 
solution that x , y , z, the kinetic energy T, and F contain no term independent 
of the arbitraries. In general, we shall avoid choosing the constant of length 
as an arbitrary so that certain functions which are only defined by their deriva- 
tives with respect to the arbitraries need no added constant independent of 
thern. 

3. Let a, a' be two constants of the cc-set, and /3, /S' two of the /3-set corres- 
ponding to angles w, w'. Denote by 8 J 8a derivatives with respect to a only in 
so far as a is contained in the coefficients, so that d/ df refers to derivatives 
when a., /3. are considered independent and 8J8f to derivatives when a., w i are 
considered independent ; evidently 8/8vJ i = 8 / 'd/3. = djdw i = d/dfi. . 
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Denote by x, y , z derivatives with respect to t. Since x — 'Z i b i dxjd^ i , we 
have 

dx dx „ dx db. 



da da l dw. da' 



(dx\ d /dx\ dx db t dx db { 

da) ~ dt\da-j l dic i da * dw. da 



d I dx 

dt 



dx _ dx db. 
5a • ». era 

in which it is to be noted that d/da, d/dt are not commutative although djdw i , 
djdt are so. 

Substituting these results in (3) and putting 



we obtain 



yr^)"" [,,/] 



__, db t db ■ 

+ *****&£ L"" W A> 

( 5 ) <2ft, 

(«,£) = [a, w] + ^i^- [w ( ,w], 

(/3, £') = [w,w']. 
Also 

(6) ^[y.*]+4[*«>1 + ^^^- - 

The quantities (f,g) must be divided into several classes, each of which 
requires separate consideration. The constants f, g may both belong to the 
a-set, or to the /3-set, or one to the a-set and the other to the /3-set. Again, they 
may both be arbitrary, or both known, or one arbitrary and one known. I shall 
only consider two known constants — one belonging to each set — so that there 
will be eight different kinds of brackets to consider, three of them being inde- 
pendent of the known constants. The extension to any number of known con- 
stants is made without difficulty. 

§ 2. Application to the Arbitrary Constants. 

4. When the constants considered in § 1 are the arbitraries of the solution, 
the right hand member of (2) is zero since F does not contain them explicitly ; 
all the (f,g) are therefore constant. But [^a, a'], [£, /3' ] consist only of 
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sums of sines owing to the forms assumed for the coordinates. Hence by § 1, 
art. 3, 

(7) (a, a') = [a, a'] = (/3, /3' ) = [w, w"\ = 0, (a, /3.) = [a, «j.], 

(8) ^ [«'.»,] -*,^ [«.»,]■ 
Also from (6), 

-[a'.u,]--^ [a, „]= — [ a ', a ] =0: 

hence 

(9) [a.,w.] = (a.,/3.) = ^.= ^^— (c^! + c 2 to s + c,w s ) (i,i=l,2,3), 

t i j 

where c 15 c 2 , c 3 are functions of the arbitrary and known constants of the a-set 
only ; they are thus completely defined in accordance with the remark at the 
end of art. 2. 

Equation (8) then gives 

db i dc t dbj dcj 

1 da da.' % da da ' 

or 

da\ l ' da J da \ l % da J 
showing that 

(10) ' Z t b t dc t =-dB, 

a perfect differential. This gives a complete definition of B if we assume that 
it contains no term independent of the arbitrary constants. 



5. Canonical systems. — The constants a. can in general be chosen arbitrarily. 
Take a. = c { . Then 

[c„w.] =(c.,£.) = l, [c,.,«).] = (c i .,^.)= (i + O- 

The choice is possible if the Jacobian of the c i with respect to the a. is not zero ; 
in the applications to be made it is easy to see that this is the case. Thus the 
c { are assumed to be independent. We then have from art. 4, 

dB 

b < = -dc.' 

where B is expressed in terms of the c. and the known constants. 

All these results are perfectly well known. The sets c., /3 { and c ; , w. are 
canonical systems, for the (f,g), \_f,g~\ ar e nothing else than Lagrange's 
brackets and the above values reduce Lagrange's equations for the variation 
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of arbitrary constants, when a disturbing function R is added to F, to the 

canonical forms 

dc { dR dftt dR 

W = d/V ~dt = ~~dc~ i '' 
and 

dc t dR' dw i dR' 

dt ~~ dtc i ' dt ~ dc i ' 

where R' = R -f B and R is supposed to be expressed in terms of the w ( , c t , t 
and known constants. 



6. Another canonical set is obtained by putting 

dc dc 2 d% 

B * = d^^ + d^/* + da;^- 

and taking a t , B t as arbitraries; this is evident from equation (9). The same 
equation shows that if we put 

dc. dc, dc. 

W * = da W i + da W * + d^ W ^ 

lit 

the a . , W form a canonical set. 

The sets c i , /3 i ; a., B ( ; a., W i have the disadvantage that when we form 
the equations of variations the time appears as a factor of the periodic terms. 
This fact prevents them from being useful for most purposes of the lunar theory 
to which application will chiefly be made. It is to be noticed nevertheless that 
when {3. or w i are taken as one set of arbitraries, the c t is the conjugate set 
which makes the whole canonical ; if a. be taken as one set then the B ( or 
W. has the same property with reference to the a. t . In general, if one set of 
constants be given the second set which renders the whole canonical is deter- 
mined uniquely ; the apparent duplicity in the above is due to the fact that W i , 
w ■ are not arbitrary constants but contain the time. 

If it be possible to choose the b i as one set of arbitraries, the conjugate sets are 

2 «.^ = 2/3 .Jl, 2.^.^=2.^.^1 ( i= l f2l 3); 

J^idb i J 3 db,' •> 3 db { ■> 3 db } 

and from (10), 

so that 

d 2 

i } 
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§ 3. Application to the Known Constants. 

7. Suppose that F contains a known constant a of the a-set and a known 
angle w i = b A t + $ t in which (the more general case) b i is a function of a and 
/3 4 belongs to the /3 set. Then applying (6) successively to the various combina- 
tions of <z , w t with the arbitraries and using (7), (9) we obtain 

W [«, «,] - sa [«\ «o] = ^ [«, «'] = 0-. 

ai [>> ™J - ^ [ «' »*] = ^ !>• "] = °' 

a r i d r i e r i 

^ K- »*] = ar K ^] -B^l a - a «] • 

4 

The first three of these give 

dC dC dC 

(li) K,™j = ^7' [ a " a o] =B_ a^' C w *' w 4] = a^: (<=i,8,8), 

and the last two then give 

U * V 4 

in which (7, , C are, owing to the forms of the square brackets as deduced from 
the expressions for the coordinates, trigonometric series with the arguments w it 
w t and with coefficients depending on a { , a ; C t consists of cosines only, G of 
sines only. Moreover, C i is only defined by its derivates with respect to a , 
« 4 , w t and may therefore contain terms independent of them, while O is only 
defined by its derivatives with respect to a { , %o i , w i and might contain terms inde- 
pendent of them. The suppositions of § 1 prevent the presence of such terms, 
so that C , C t are fully defined. 

If the solution of the original equations has been obtained algebraically we 
have then sufficient material to find C 4 , C . It may happen, however, that 
numerical values have been given to the constants in forming the solution, and 
for many purposes it is useful to possess equations for C 4 , C which will enable 
us to find their values in this case. It will be shown that they satisfy ordinary 
differential equations of the first order. 



340 E. W. BROWN : ON THE VARIATION [July 

These equations may be obtained directly from the results of arts. 3, 4, 7 
without using the original equations of motion. But as they involve d'F/da , 
d'F/dw t , it is much easier to obtain them directly from the results of arts. 4, 7 
and the equations of motion. This latter method, which also furnishes the vis 
viva equation, will be followed. 

In this section, unless explicitly stated otherwise, the suffixes i,j receive only 
the values 1, 2, 3, that is, the letters to which they are attached are arbitraries 
of the solution. 



8. The Vis Viva equation. Since 

r t dc J r i dC i 

dQ_ dQ dQ 
dt ~ t0i dv>~ t + *dto t ' 

where Q contains the time only through its presence in the angles w t , w 4 , we 
have 

dc- dC a (dx . ,.dx\ d ,m et\ 

2.6. -j± + 6, -j-* = S[ 5— x — x ^r— ) = ^- ( T — F) . 

' } da^ * da { \ da t da t ) da t K ' 

But 2 . 6 . dc . = — dJB . Therefore, since there are no terms independent of the 
arbitrary constants and 6 4 is a known constant, 

(13) T-F-b 4 C 4 =-B. 

9. The equation for C 4 . Differentiating (13) with respect to the time, 

b ^dt-Jt( T -^- S ^-^- h ^ 
or 

(14) d_C i= _d'F^ 
v ' dt dw, 



The variable C 4 is therefore similar to the L of M. Poincae^ in his paper 
Sur les equations de la lune. * We might have defined C 4 by means of this 
equation and have then reversed the processes of this article and obtained the 
first and third of equations (11). As 4 can contain a constant term, the defi- 
nition is not complete, but will be made so by the equation obtained in art. 9. 



Bulletin Astronomique, vol. 17 (1900), p. 167. 
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10. The equation for C Q . — Multiply the first of equations (12) by b t , the 
second by 6 4 and form the sum. We obtain 

2,6, 



da J da 


dC /Sac -, dx \ 
at \ oa n oa ) 




DT DF d'F 




" da ~ da da 




5 dB d'F 

-da^)-d^ + l3or 



by (13). In this equation B is supposed to be expressed in terms of a t , a^. 
Let dBJda denote the derivative of B with respect to a when B is expressed 
in terms of c 15 c 2 , c 3 , a , so that 

dB dB e dB o dc l dB c dc t 

da da Q l dc t da da l * da Q ' 

Substituting in the former equation and rearranging, we obtain 

nn dC db 4 _dB^_ d'F 

{ } ~dt' + 4 ffc„ da ~ 'B\ ' 

the equation required. 

Since C contains no constant term, (15) gives it completely when C 4 has 
been obtained. We already know C 4 from (14) except its constant term. But 
since C contains no term proportional to the time, the constant term in O i is 
given by 

where ( • ■ )° denotes the constant term in the expansion of the enclosed function 
in trigonometric series. 

If b t does not depend on a , it is easily seen that the equation has its second 
term absent and C can be found independently of C i . 

But now the constant part of C 4 cannot be determined by (15)°. The diffi- 
culty is turned by following the investigation of this article with b 4 as the con- 
stant instead of a . Then as b t occurs explicitly in i^only through its presence 
in w i the equation (15) reduces to 

giving 

(OJ '"■ 



db i 



The C of this equation is, of course, different from that of (15). 

Trans. Am. Math. Soc. 33 
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There may be other known constants a present in F. If b i depends on a 
then there will be an equation for C corresponding to (15). If b t is indepen- 
dent of a , the equation satisfied by the corresponding C j is 

(16) ^_^ = _?Z 

^ 10 > dt da Ba ' 



where C consists of sines only. It will be convenient in what follows to sup- 
pose that any other constants are all of the latter type and therefore may either 
be treated as numerical constants in which case the corresponding equation (16) 
will be unnecessary, or as algebraic constants in which case it may be necessary 
to use the equation. Both cases may arise in the method of the variation of 
known constants which is developed below. 

11. Suppose that F is homogeneous of order q with respect to x, y, z, and 
a certain constant A present in F, there being also present an absolute constant 
(e. g., the " gravitation " constant) which restores the proper physical dimensions 
of the equations of motion. 

Then we have, symbolically, 

F, B, T=q (length). 
The equations of motion show that 

time = (1 — iq) (length), 
and since [c i ,w < ]=l,6= — 1 (time), 



c, 



(1 + iq) (length) = ^(ft,), F, B, T= ^(c,). 



The choice of a has been quite arbitrary. When it is a function of b t take 
it so that it is of the same dimensions as the c. and further such that 

<K _ q + 2 a 

db^~ q — 2 b' t ' 

As A is merely denned as a length, choose it so that 

dA 2 da 

~A~ = q + 2 ~a • 

Thus A, a are defined in terms of 6 4 , but the manner in which they occur in 
.Pis still arbitrary except with respect to their dimensions. 

Since B will be homogeneous with respect to the c., a we have 

dB c 2q dB 2q 

a„ ~7^ = — S> S — z.c, -j— = — -"o B + 2,. i b i c . 
da t) q+2 l dc t q + 2 ' l 
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12. Multiply (15) by q + 2 and apply the results for dA , da in art. 11. 
We obtain 

where F is supposed to be expressed in terms of A instead of a , 

Now F is homogeneous of order q with respect to x , y , z , A . Hence, mul- 
tiplying equations (1) by 2x , 2y , 2z , and adding their sum to the last equation, 
we obtain 

^-4T-(q+2)a ^-(q-2)b i C i +( q +2)a -^ = 2qF, 

which, by (15) and the last result in art. 11, becomes 

r n fin 

^r*-(q+2)a ^-(q-2)(T-F+B) + 2 q B 

+ (q + 2)?Lb i c i =2qF+AT, 
or, reducing, 

(17) T+F-B-b^-b^-b.c^j^^-^-a^. 

Denote by zero index the constant terms in the expansion of these functions. 
Then, since the right-hand member contains no constant term, 

(18) T° + F° = B + b lCl + b 2 c 2 + b^, 
a result independent of q . 

13. The equation (18) gives a useful formula for finding the c { when the 
solution has been obtained with any set of arbitraries and even when numerical 
values have been given to the a-set. For 

c i= -^(B + Zb i c i ) = ~(f» + F<>) 

= S [_*W i + X ^\ = S l X dv i + dt{ X 8b i )\ < 
or 

[dx . Dy . dz 1° 

Define a new constant c 4 by a similar equation for i = 4 . We then have 

6 1 c 1 + 5 2 c 2 + 6 3 c 3 +5 4 c 4 = 2 7 7 % 
or by (18), 

T°-F°-b i c i =-B. 

4 4 

Comparing with the vis viva equation, we obtain 

(20) c 4 = {CJ. 
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§ 4. Canonical Equations of Motion. 
14. The equations (14), (15) may be written 



(21) *-& — **, <(o.-, d *>)-^™°£ 

v ' dt otv i dt\ ° oa J oa g * da t 



' dx d(f> 


dx 


d(f> 




~dt = dx ' 


dt = 


~dx ; " 


' f ' ' ' 


dC 4 d<f> 


dw i 


d<j> 




dt ~ dw. ' 

^ 4 


~dt = 


~dCj 




da 


B<f> 


dW 


d(f> 


~di = 


dW' 


dt 


~da n 



These with the equations of motion (1) suggest the following transformation. 
Put 

2T=x* + y* + z\ W =-C -b g t + f3 , = T - F — 6 4 C 4 , 

where /3 is a constant. If we consider x , y , x , (7, , a and x , y, z, iv 4 , W as 
independent variables and <f> expressed as a function of them, the equations (1) 
and (21) may be written in the canonical form 



(22) 



(23) 

for <f> does not contain W explicitly, « is a constant, and dw i jdt= b±. The 
same equations hold if b t is not a function of a g . 

We may either use the first four or all five pairs of equations and the solu- 
tions of such canonical systems are well known. The last pair gives a g = 
constant and defines W . It will be more convenient, however, to treat « as a 
given constant and use only the first four equations in the following work. 

When (22) are solved canonieally in terms of the constants c. and the angles 
w J (i=l, 2, 3, 4) we find 

TV 7 dR 

<p = constant = — B , o. = 5— . 

^ J dc. 

j 

The previous articles show that B' = B + b i c A where B is independent of c 4 
[the assumption that the" constant parts of (f>, etc., contain no parts independent 
of c n c 2 , c 3 no longer holds]. Hence, as before, b t = — dB i /oc i . Also the 
expressions for the coSrdinates are independent of c t . 

15. All the preceding investigations might have been developed by starting 
with these canonical equations and solving them by means of Jacobi's method, 
obtaining the results in the inverse order. I have preferred to follow the more 
natural if somewhat longer development as it seems to make the deduction of 
some of the results rather more clear in view of the applications which follow. 
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If there be angles other than w i and constants other than « present in F, 
similar results may be obtained. In fact, the extension of this method to any 
number of variables, angles and constants will be sufficiently obvious. To each 
angle and constant corresponds respectively a pair of canonical equations ; to <f> 
must be added one term for each angle present. When this has been done, <j>, 
considered as a function of all the variables, is independent of t explicitly and 
$ = const, is an integral of the equations. 

§ 5. Variation of Constants. 

16. Variation of arbitrary constants. — Let a new function ty be added to 
(j), then, using the canonical set c,, c 2 , c 3 , c 4 , w l , w a , w 3 , W 4 , the equations sat- 
isfied by them are 

where ^ and 5' are expressed in terms of c. , w y , of the given constants and the 
time, B' being a function of the c . and the constants present in </> . 

The method appears to fail if yjr contains none of the arbitrages c , tv.. 
But the equations of motion depend on the derivatives of \]r with respect to the 
variables ; we can therefore assume that yjr is not independent of the arbitraries. 

If i|r is independent of the w } , the c 4 are constant and the equations for w } 
are easily solved. If ijr contains a small factor a and a term with argument 
bt + ft, where b is small and /3 independent of the /3 4 , the corresponding terms 
in w t will contain divisors b, ¥ja, b 3 /a 2 ,' •••; when /3 is not independent of 
/3. the first divisor will in general be o 2 . 

17. Variation of the given constants. Let -^ be due to a small addition Sa 
to a , one of the constants present in F where 8a Q may be a constant or a func- 
tion of the time. Then using the canonical equations (22) and assuming that 
Sa produces variations 6c, Bio. and that squares of the variations may be 
neglected, 

ylr = 2. ( V 1 ec. + a — dw, + -_— 6a. 1 = — 2, — ?— 6c. -j— oa n , 

r > \ dc. 3 dw. 3 da a J ■> dc j ' da Q °' 

since 2?' is independent of the w. . 
The equations of variations are 

d . df dB' d dB' d 

-r- oc . = — ~- = 2 ., -r: — _ — oc,, + -^ — ^— oa. , 
di! ' dw. 3 ac, cm 3 da v ow j ° 

d , * n a ^ j. d (^ dB \ , dB \ \ , A 

(i,/ = l,2, 3,4). 
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As squares of the variations are neglected we can substitute constant values 
for the coefficients of the variations in the right hand members and therefore 

-77OC = — 2.,0.,^r dC, + ~t ■ = 02 n , 

dt > J > dw. ■> ^ da n dw. °' 

3 v 3 



d 

dt 



(w, + Szo,) = ^ (?..,b.,8c.,--^ 8a o y b.. 



18. Suppose that Ba is independent of the w., c,. The last terms of the 
equations for Sc . then vanish and they are satisfied by 

8c. = 0. 

3 

Substituting in the second set we find, since B' = B + b i c i , 
d r , s \ db J s , j. 

-y- ( W. + OW.) = -T+- O0T.+ 0.. 

eft v •> ^ ' ' da 0^-3 

That is, in order to obtain the new value of w , all that is necessary is to 
substitute a + ba a for a in the expression for b, expressed as a function of 

C l» C 2' C 3' %• 

Suppose that the original solution had been expressed in terms of any arbi- 
traries a, , a 2 , a 3 of the a-set, and of the known constant a ; the c t having been 
obtained in terms of them by the method of art. 13. Then the variations of <x t 
are found by solving the equations 

dc t dc h dc t dc t 

da, **+ da^ + d-a.^ = ~ da*** (l = 1 ' 2 ' 3) ' 

for 8a,, ba 2 , 8a s . The changes in the coefficients of the periodic terms are 
obtained by substituting these changes of the constants therein. But the 
changes in the b t are given by 

., db. . db, . db. . db. . 

8b, = ^ 8a, + ^ 8a 2 + ^ Ba, + ^ S« (« = i, ., 3), 

that is to say, we must also substitute the change of a as well as those of 
a,, a x , a 3 in the b i . 

The same results hold if b. is not a function of a. . 

19. If ty be due to a variation bw i of «> 4 , we have in a similar manner, 

■*ir = 2 , b . Sc . ; 

~ j 3 j ' 

and the equations of variations are 
Jt 8c. = - 2^.,~ &y, !( W , + fa,- ) - - |- (S-.fU,,) + 6,. 



1903] OF CONSTANTS IN DYNAMICAL EQUATIONS 347 

Suppose that the variation of w i is given. The left hand member of the last 
of the second set, namely, 

-j- Sw, = 2 ., 6 ., v— 8c., 
dt 4 ■> 3 dc, •> 

is known. 

The value of 8w i may consist of periodic and non-periodic terms : suppose 
it to be independent of the angles w 1 , w 2 , tv 3 and of c l , c 2 , c 3 . Let us try 
if Scj = 0, Sc 2 = 0, Sc 3 = and Sc 4 independent of w 1 , w 2 , w. i can satisfy the 
equations. The first three of the first pair are satisfied and the others become 

jt^ = - b ^L^ dt Sw * =h *dc 4 Sc » jt Sw < =b ^ Sc * < ! ' =i ' 2 ' 3) - 

4 4 * 

If Sw 4 = qt + <?> where q, q are given constants, the equations are satisfied 
by Bc i = qc i /b i , 8w i = 0. It does not appear to be possible in general to 
satisfy the equations with other values for $w i . 

§ 6. Applications to the Lunar Theoey. 

20. The results contained in the previous sections are more general in their 
character than is necessary for the applications to most problems in celestial 
mechanics and in particular for that one which is known as the lunar theory. 
Without going into details it is sufficient to state here that the few restrictions 
imposed are satisfied in these applications. 

Let x, y, z, r denote the coordinates and distance of the moon relatively to 
the earth's center ; x', y\ »', r\ those of the sun relatively to the center of 
mass of the earth and moon ; \x = E + M, the sum of the masses of the earth 
and moon ; m the mass of the sun ; let the three bodies be considered to attract 
in the same way as particles of the same masses situated at their centers of 
mass respectively. Also suppose that the motion of the sun is known and that 
all other influences are neglected. Then it is well known that 

F= ~ + m! \r — 2r(xx' + yy' + zz') +r' } _ * — m'(xx'-\- yy' + zz')r' , 

where, after the substitution of polar coordinates for x', y', z and expansion in 
powers of 1/r', the term 1/r' is dropped and the coefficient of {r')~ p is multi- 
plied by 

_ / e y- 1 / - m y- 1 
K p - \e~+~m) ~ \e + m) ■ 

We can therefore put 

E=^ + m'±,c p l}' p (r')-r-\ 

where R' p is a homogeneous function of order ^3 with respect to x, y, z, and also 
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contains trigonometric functions of the angular position of the sun. Further, if 
n , a be the mean motion and mean distance of the sun, m + /x = n a , or, 
putting X = /x/(m' + /jl), 

m = n a' (1 — X). 

We shall now suppose that the plane of the sun's orbit is fixed and that its 
motion can always be expressed in terms of n , e' and the angle n't + /3 t ; when 
its motion is elliptic n' , /3 4 , e' are constants. Properly, there are two angles 
n't + e' , n't + e' — tr' , but well known results show that we may treat them as 
though 7r' were zero and they also indicate where ir' has to be inserted in the 
final results. Thus we can put 

R r 



p=2 CI 



r 



where 7? ) is homogeneous and of order p relatively to x, y, z and contains k , 
e' and trigonometric functions of n't + e' . 

21. We have from the equations of motion, 

I C ^- - %T= 8x~ = -^ + n'\l - X^pR^a')-^ { P =2, 3, •••)• 
Eliminating T between this equation and (17) we find, since q = — 1 here, 

(24) | £ -$n'\l-\)Z p (p -2 ) R v {a' )-*+*- B - Z t b t c t = | ^ - a -J' . 
The right hand member contains no constant term. Hence 

(25) J e + 5 1 c 1 + 5 2 c 2 + 6 8 c 3 = |^y-^' 2 (l-X)2 ) (F-2)(i? ji ) u («')- +2 . 

The equations (14), (15) show that dC /dt is a quantity whose lowest order 
is n /a u . As it is usual in the lunar theory to expand in powers of n and \Ja , 
it will be more convenient to put a g C — C' so that dC' /dt is a quantity whose 
lowest order is n" and is expansible like the coordinates. We therefore have 

d C d 2 u. 

(26) - d f = ^ 2 ^)-^+l-n'\l-X)2X P -2)R p (a')-^+f 

(P = 3,4,---), 

f being a constant so chosen as to make the constant part of the right-hand 
member zero. 

22. The properties of equations (25), (26) which are remarkable and valuable 
for practical purposes arise principally in those investigations which permit us 
to neglect the ratio of the parallaxes of the sun and moon. This ratio is 
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about j^ and for most of the purposes which require the value of C' u , i. e., 
in the calculation of the indirect action of a planet (art. 24 below), it is possible 
to neglect it. But of the R , the first term i? 2 is the only one which is inde- 
pendent of it and the factor p — 2 causes R 2 to disappear in both the equations. 
The equation (26) then reduces to 

Another observation may be made. If the object in view is such that we 
only require those terms in C' Q which are of long period relatively to the month, 
i. e., those in which the coefficients of t in the arguments are small, the first 
term of the right-hand member will give rise to coefficients which are very small 
compared with those arising from the second member ; in such cases we can 
generally put 

---•(?)♦•(?)'• 



dC 



dt 
without committing sensible errors 



'o 



23. In nearly all cases of the indirect action of a planet, we can neglect the 
square of the ratio of the parallaxes. But it is well known that the constant 
parts of the functions we have been dealing with are expansible in powers of the 
square of this ratio. Hence neglecting the square of the ratio of the parallaxes, 



B + b.c. + Jc, + b,c. 



»(?)'■ 



This is the result which I obtained in an earlier paper * and which I used to 
prove Adams' well known theorems and to obtain the secular accelerations 
directly from the expression for the parallax of the moon. It is true that the 
square of the ratio of the parallaxes was neglected, but it was easy to see from 
the few terms which Delaunay gave in a literal form that they could exert no 
sensible effect on the result. If doubts might still arise, owing to slow conver- 
gence, it is not difficult to remove them by means of the numerical results of my 
lunar theory taken in conjunction with equation (25). 

The methods used here differ in form from those which were adopted in the 
paper just referred to. The essential feature of the earlier method was really 
that of Delaunay, that is to say, the elliptic orbit was used as a first approxima- 
tion and the variation of arbitrary constants applied to it somewhat in Delau- 
nay's manner. Here we are theoretically quite independent of the elliptic orbit, 
and any set of arbitrary constants which possesses the property of enabling 

*See footnote, p. 334 above. 
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us to expand in trigonometric series may be used. This has the practical advan- 
tage of permitting the calculation of c,, c 2 , c 3 , by means of art. 13 from the 
results of any theory in whatever form it may be expressed and whatever be the 
arbitrary constants used. Provisionally it may be best to except the theory of 
Hansen, although it is probable that the expression xBx + i/By + zSz, which 
is the principal part of equation (19), might be transformed so as to be used in 
this case also Without serious difficulty. 

24. The results of arts. 15, 16 include Newcomb's theorem concerning the 
secular accelerations and the indirect action of a planet * as a particular case. 
In this valuable memoir he considered the action of a planet as a case of the 
disturbed motion of three bodies, that is to say the problem of the rela- 
tive motion of the earth, moon and sun is considered to have been solved and 
the question was then to find out what variable values the 12 arbitrary con- 
stants which arise in the solution must have in order that the solution may 
retain the same form when another disturbing cause acts. He showed that if 
we first consider the orbit of the sun as elliptic and if the action of a planet 
causes variations Be', Sn' in e , n , then the new non-periodic parts of n, tt 1 , 6 1 
(the mean motions of the moon, its perigee and node) are found by solving the 
equations 

-^~~ on + -=- be + -=- 07 -f ■=— ,* bn + ^.—7 be = ( * = 1 , 2 , 3 ) 
on de oj en oe 

for Sn, Be, By, regarding Sn', Be' as known, and substituting the variable values 
of 11, e, 7, n , e in the expressions of n, 77^, 8 X in terms of these five quantities. 
He remarks, " This curious theorem may embody some principle applicable to 
the disturbed motion of three bodies, which has not yet been fully mastered." 

It appears from the investigation of this paper that the theorem is independent 
of the particular form in which the solar coordinates enter into the lunar dis- 
turbing function and further, that it is true for any of the constants which may 
be present in this disturbing function, provided always that the given variations 
of the known constants are independent of the lunar elements. Finally, it is 
applicable to any dynamical problem in which the coordinates satisfy the speci- 
fied conditions and to any constants which may enter into the force function in 
the manner specified at the outset. 

Havkeford College, 
April 20, 1903. 

* Action of the Planets on the 3f nt»i, Washington Astronomical Papers, vol. 5 (189-1), 
p. 191. 



